Replacement of the Gamma by the Delta variant in Brazil: impact of lineage displacement on the ongoing pandemic (preprint)

The COVID-19 epidemic in Brazil was driven mainly by the spread of Gamma (P.1), a locally emerged Variant of Concern (VOC) that was first detected in early January 2021. This variant was estimated to be responsible for more than 96% of cases reported between January and June 2021, being associated with increased transmissibility and disease severity, a reduction in neutralization antibodies and effectiveness of treatments or vaccines, as well as diagnostic detection failure. Here we show that, following several importations predominantly from the USA, the Delta variant rapidly replaced Gamma after July 2021. However, in contrast to what was seen in other countries, the rapid spread of Delta did not lead to a large increase in the number of cases and deaths reported in Brazil. We suggest that this was likely due to the relatively successful early vaccination campaign coupled with natural immunity acquired following prior infection with Gamma. Our data reinforces reports of the increased transmissibility of the Delta variant and, considering the increasing concern due to the recently identified Omicron variant, argues for the necessity to strengthen genomic monitoring on a national level to quickly detect and curb the emergence and spread of other VOCs that might threaten global health.

Genomic epidemiology reveals how restriction measures shaped the SARS-CoV-2 epidemic in Brazil (preprint)

Brazil has experienced some of the highest numbers of COVID-19 infections and deaths globally and made Latin America a pandemic epicenter from May 2021. Although SARS-CoV-2 established sustained transmission in Brazil early in the pandemic, important gaps remain in our understanding of local virus transmission dynamics. Here, we describe the genomic epidemiology of SARS-CoV-2 using near-full genomes sampled from 27 Brazilian states and an adjacent country - Paraguay. We show that the early stage of the pandemic in Brazil was characterised by the co-circulation of multiple viral lineages, linked to multiple importations predominantly from Europe, and subsequently characterized by large local transmission clusters. As the epidemic progressed, the absence of effective restriction measures led to the local emergence and international spread of Variants of Concern (VOC) and under monitoring (VUM), including the Gamma (P.1) and Zeta (P.2) variants. In addition, we provide a preliminary genomic overview of the epidemic in Paraguay, showing evidence of importation from Brazil. These data reinforce the need for the implementation of widespread genomic surveillance in South America as a toolkit for pandemic monitoring and providing a means to follow the real-time spread of emerging SARS-CoV-2 variants with possible implications for public health and immunization strategies.

Asymptomatic SARS-CoV-2 infection and the demography of COVID-19 (preprint)

Asymptomatic individuals carrying SARS-CoV-2 can transmit the virus and contribute to outbreaks of COVID-19, but it is not yet clear how the proportion of asymptomatic infections varies by age and geographic location. Here we use detailed surveillance data gathered during COVID-19 resurgences in six cities of China at the beginning of 2021 to investigate this question. Data were collected by multiple rounds of city-wide PCR test with detailed contact tracing, where each patient was monitored for symptoms through the whole course of infection. We find that the proportion of asymptomatic infections declines with age (coefficient =-0.006, P<0.01), falling from 56% in age group 0-9 years to 12% in age group >60 years. Using an age-stratified compartment model, we show that this age-dependent asymptomatic pattern together with the age distribution of overall cases can explain most of the geographic differences in reported asymptomatic proportions. Combined with demography and contact matrices from other countries worldwide, we estimate that a maximum of 22%-55% of SARS-CoV-2 infections would come from asymptomatic cases in an uncontrolled epidemic based on asymptomatic proportions in China. Our analysis suggests that flare-ups of COVID-19 are likely if only adults are vaccinated and that surveillance and possibly control measures among children will be still needed in the future to contain epidemic resurgence.

Fatal COVID-19 outcomes are associated with an antibody response targeting epitopes shared with endemic coronaviruses (preprint)

It is unclear whether prior endemic coronavirus infections affect COVID-19 severity. Here, we show that in cases of fatal COVID-19, antibody responses to the SARS-COV-2 spike are directed against epitopes shared with endemic beta-coronaviruses in the S2 subunit of the SARS-CoV-2 spike protein. This immune response is associated with the compromised production of a de novo SARS-CoV-2 spike response among individuals with fatal COVID-19 outcomes.

Analysis and visualization of epidemics on the timescale of burden: derivation and application of Epidemic Resistance Lines (ERLs) to COVID-19 outbreaks in the US (preprint)

The 2020 COVID-19 pandemic produced thousands of well-quantified epidemics in counties, states, and countries around the world. Comparing the dynamics and outcomes of these nested epidemics could improve our understanding of the efficacy of non-pharmaceutical interventions (NPIs) and help managers with risk assessment across multiple geographic levels. However, cross-outbreak comparisons are challenging due to their variable dates of introduction of the SARS-CoV-2 virus, rates of transmission, case detection rates, and asynchronous and diverse management interventions. Here, we present a graphical method for comparing ongoing COVID-19 epidemics by using disease burden as a natural timescale for comparison. Trajectories of growth rates of cases over the timescale of lagged deaths per-capita produces coherent visual comparisons of epidemics that are otherwise incoherent and asynchronous in the timescale of calendar dates or incomparable using non-stationary measures of burden such as cases. Applied to US COVID-19 outbreaks at the county and state level, this approach reveals lockdowns reducing transmission at fewer deaths per-capita early in the epidemic, reopenings causing resurgent summer epidemics, and peaks that while separated in time and place actually occur at points of similar per-capita deaths. Our method uses early and minimally mitigated epidemics, like that in NYC in March-April 2020 and Sweden in later 2020, to define what we call “epidemic resistance lines” (ERLs) or hypothesized upper bounds of epidemic speed and burden. ERLs from less-mitigated epidemics allow benchmarking of resurgent summer epidemics in the US. In particular, the unmitigated NYC epidemic resistance line appears to bound the growth rates of 3,000 US counties and funnel growth rates across counties to their peaks where growth rates equal zero in the fall and winter of 2020. Corroboration of upper-bounds on epidemic trajectories allowed early predictions of mortality burden for unmitigated COVID-19 epidemics in these populations, predictions that were more accurate for counties in states without mask-wearing mandates. We discuss how this method could be used for future epidemics, including seasonal epidemics caused by influenza or ongoing epidemics caused by new SARS-CoV-2 variants. Press Summary Why, despite no statewide mask-wearing mandates or other restrictions like restaurant closures, did South Dakota’s COVID-19 epidemic peak not in January, when seasonal forcing wanes, but in early November? Why are we not seeing a resurgent epidemic in Florida or Texas, where non-pharmaceutical interventions have been relaxed for months? How can we compare the current outbreak in India with other countries’ epidemics to contextualize the speed of the Indian outbreak and estimate the potential loss of life? We have developed a new method of visualizing epidemics in progress that can help to compare distinct COVID-19 outbreaks to understand, in specific cases like South Dakota, why they peaked when they did. The “when” in this case does not refer to prediction of a calendar date, but rather a point in the accumulation of deaths in a given locale due to the disease in question. The method presented in this paper therefore essentially uses population-based burden of disease as a timescale for measuring epidemics. Just as the age of a car can be measured in years or miles, the age of a COVID-19 epidemic can be measured in days or deaths per-capita. Plotting growth rates of cases as a function of per-capita deaths 11 days later produces a real-time visual comparison of epidemics that are otherwise asynchronous in time. This approach permits both direct comparison across local outbreaks that may be disparate in time and/or place, as well as benchmarking of any outbreak against known exemplars of archetypal response strategies, such as New York City’s unmitigated urban outbreak in Spring 2020 and Sweden’s uncontained summer 2020 epidemic. Whether comparing the speed of resurgent outbreaks following relaxation in US states like Florida or the peak mortality burden in fall outbreaks across thousands of US counties with and without statewide mask-wearing mandates, this method offers a simple, intuitive tool for real-time monitoring and prediction capability connecting epidemic speed, burden, and management interventions. While our findings point to compelling epidemiological hypotheses for peaks in less-regulated states, future work is needed to confirm and extend our results predicting mortality burden at the peak of confirmed cases in the ongoing and evolving COVID-19 pandemic.

Early genomic detection of SARS-CoV-2 P.1 variant in Northeast Brazil (preprint)

Tracking the spread of SARS-CoV-2 variants of concern is crucial to inform public health efforts and control the ongoing pandemic. Here, we report genetic evidence for circulation of the P.1 variant in Northeast Brazil. We advocate for increased active surveillance to ensure adequate control of this variant throughout the country. Article Summary Line Active genomic surveillance of SARS- CoV-2 suspected cases from recent travelers reveals the circulation of the P1 variant of concern in Bahia state, Northeast Brazil.

Emergence and rapid spread of a new severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) lineage with multiple spike mutations in South Africa (preprint)

Continued uncontrolled transmission of the severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) in many parts of the world is creating the conditions for significant virus evolution. Here, we describe a new SARS-CoV-2 lineage (501Y.V2) characterised by eight lineage-defining mutations in the spike protein, including three at important residues in the receptor-binding domain (K417N, E484K and N501Y) that may have functional significance. This lineage emerged in South Africa after the first epidemic wave in a severely affected metropolitan area, Nelson Mandela Bay, located on the coast of the Eastern Cape Province. This lineage spread rapidly, becoming within weeks the dominant lineage in the Eastern Cape and Western Cape Provinces. Whilst the full significance of the mutations is yet to be determined, the genomic data, showing the rapid displacement of other lineages, suggest that this lineage may be associated with increased transmissibility.

SARS-CoV-2 shifting transmission dynamics and hidden reservoirs limited the efficacy of public health interventions in Italy (preprint)

We investigated SARS-CoV-2 transmission dynamics in Italy, one of the countries hit hardest by the pandemic, using phylodynamic analysis of viral genetic and epidemiological data. We observed the co-circulation of at least 13 different SARS-CoV-2 lineages over time, which were linked to multiple importations and characterized by large transmission clusters concomitant with a high number of infections. Subsequent implementation of a three-phase nationwide lockdown strategy greatly reduced infection numbers and hospitalizations. Yet we present evidence of sustained viral spread among sporadic clusters acting as "hidden reservoirs" during summer 2020. Mathematical modelling shows that increased mobility among residents eventually catalyzed the coalescence of such clusters, thus driving up the number of infections and initiating a new epidemic wave. Our results suggest that the efficacy of public health interventions is, ultimately, limited by the size and structure of epidemic reservoirs, which may warrant prioritization during vaccine deployment.

Major new lineages of SARS-CoV-2 emerge and spread in South Africa during lockdown. (preprint)

In March 2020, the first cases of COVID-19 were reported in South Africa. The epidemic spread very fast despite an early and extreme lockdown and infected over 600,000 people, by far the highest number of infections in an African country. To rapidly understand the spread of SARS-CoV-2 in South Africa, we formed the Network for Genomics Surveillance in South Africa (NGS-SA). Here, we analyze 1,365 high quality whole genomes and identify 16 new lineages of SARS-CoV-2. Most of these unique lineages have mutations that are found hardly anywhere else in the world. We also show that three lineages spread widely in South Africa and contributed to ~42% of all of the infections in the country. This included the first identified C lineage of SARS-CoV-2, C.1, which has 16 mutations as compared with the original Wuhan sequence. C.1 was the most geographically widespread lineage in South Africa, causing infections in multiple provinces and in all of the eleven districts in KwaZulu-Natal (KZN), the most sampled province. Interestingly, the first South-African specific lineage, B.1.106, which was identified in April 2020, became extinct after nosocomial outbreaks were controlled. Our findings show that genomic surveillance can be implemented on a large scale in Africa to identify and control the spread of SARS-CoV-2.

Potential impact of individual exposure histories to endemic human coronaviruses on age-dependence in severity of COVID-19 (preprint)

Cross-reactivity to SARS-CoV-2 from previous exposure to endemic coronaviruses (eHCoV) is gaining increasing attention as a possible driver of both protection against infection and severity of COVID-19 disease. Here, we use a stochastic individual-based model to show that heterogeneities in individual exposure histories to endemic coronaviruses are able to explain observed age patterns of hospitalisation due to COVID-19 in EU/EEA countries and the UK, provided there is (i) a decrease in cross-protection to SARS-CoV-2 with the number of eHCoV exposures and (ii) an increase in potential disease severity with number of eHCoV exposures or as a result of immune senescence. We also show that variation in health care capacity and testing efforts is compatible with country-specific differences in hospitalisation rates. Our findings call for further research on the role of cross-reactivity to endemic coronaviruses and highlight potential challenges arising from heterogeneous health care capacity and testing.

The impact of host resistance on cumulative mortality and the threshold of herd immunity for SARS-CoV-2 (preprint)

It is widely believed that the herd immunity threshold (HIT) required to prevent a resurgence of SARS-CoV-2 is in excess of 50% for any epidemiological setting. Here, we demonstrate that HIT may be greatly reduced if a fraction of the population is unable to transmit the virus due to innate resistance or cross-protection from exposure to seasonal coronaviruses. The drop in HIT is proportional to the fraction of the population resistant only when that fraction is effectively segregated from the general population; however, when mixing is random, the drop in HIT is more precipitous. Significant reductions in expected mortality can also be observed in settings where a fraction of the population is resistant to infection. These results help to explain the large degree of regional variation observed in seroprevalence and cumulative deaths and suggest that sufficient herd-immunity may already be in place to substantially mitigate a potential second wave.

Early transmission of SARS-CoV-2 in South Africa: An epidemiological and phylogenetic report (preprint)

Background: The emergence of a novel coronavirus, SARS-CoV-2, in December 2019, progressed to become a world pandemic in a few months and reached South Africa at the beginning of March. To investigate introduction and understand the early transmission dynamics of the virus, we formed the South African Network for Genomics Surveillance of COVID (SANGS_COVID), a network of ten government and university laboratories. Here, we present the first results of this effort, which is a molecular epidemiological study of the first twenty-one SARS-CoV-2 whole genomes sampled in the first port of entry, KwaZulu-Natal (KZN), during the first month of the epidemic. By combining this with calculations of the effective reproduction number (R), we aim to shed light on the patterns of infections that define the epidemic in South Africa. Methods: R was calculated using positive cases and deaths from reports provided by the four major provinces. Molecular epidemiology investigation involved sequencing viral genomes from patients in KZN using ARCTIC protocols and assembling whole genomes using meticulous alignment methods. Phylogenetic analysis was performed using maximum likelihood (ML) and Bayesian trees, lineage classification and molecular clock calculations. Findings: The epidemic in South Africa has been very heterogeneous. Two of the largest provinces, Gauteng, home of the two large metropolis Johannesburg and Pretoria, and KwaZulu-Natal, home of the third largest city in the country Durban, had a slow growth rate on the number of detected cases. Whereas, Western Cape, home of Cape Town, and the Eastern Cape provinces the epidemic is spreading fast. Our estimates of transmission potential for South Africa suggest a decreasing transmission potential towards R=1 since the first cases and deaths have been reported. However, between 06 May and 18 May 2020, we estimate that R was on average 1.39 (1.04 - 2.15, 95% CI). We also demonstrate that early transmission in KZN, and most probably in all main regions of SA, was associated with multiple international introductions and dominated by lineages B1 and B. The study also provides evidence for locally acquired infections in a hospital in Durban within the first month of the epidemic, which inflated early mortality in KZN. Interpretation: This first report of SANGS_COVID consortium focuses on understanding the epidemic heterogeneity and introduction of SARS-CoV-2 strains in the first month of the epidemic in South Africa. The early introduction of SARS-CoV-2 in KZN included caused a localized outbreak in a hospital, provides potential explanations for the initially high death rates in the province. The current high rate of transmission of COVID-19 in the Western Cape and Eastern Cape highlights the crucial need to strength local genomic surveillance in South Africa.

The ongoing COVID-19 epidemic in Minas Gerais, Brazil: insights from epidemiological data and SARS-CoV-2 whole genome sequencing. (preprint)

The recent emergence of a previously unknown coronavirus (SARS-CoV-2), first confirmed in the city of Wuhan in China in December 2019, has caused serious public health and economic issues due to its rapid dissemination worldwide. Although 61,888 confirmed cases had been reported in Brazil by 28 April 2020, little was known about the SARS-CoV-2 epidemic in the country. To better understand the recent epidemic in the second most populous state in southeast Brazil (Minas Gerais, MG), we looked at existing epidemiological data from 3 states and sequenced 40 complete genomes from MG cases using Nanopore. We found evidence of multiple independent introductions from outside MG, both from genome analyses and the overly dispersed distribution of reported cases and deaths. Epidemiological estimates of the reproductive number using different data sources and theoretical assumptions all suggest a reduction in transmission potential since the first reported case, but potential for sustained transmission in the near future. The estimated date of introduction in Brazil was consistent with epidemiological data from the first case of a returning-traveler from Lombardia, Italy. These findings highlight the unique reality of MGs epidemic and reinforce the need for real-time and continued genomic surveillance strategies as a way of understanding and therefore preparing against the epidemic spread of emerging viral pathogens.

Neutralising antibodies to SARS coronavirus 2 in Scottish blood donors - a pilot study of the value of serology to determine population exposure (preprint)

BackgroundThe progression and geographical distribution of SARS coronavirus 2 (SARS-CoV-2) infection in the UK and elsewhere is unknown because typically only symptomatic individuals are diagnosed. We performed a serological study of blood donors in Scotland between the 17th of March and the 18th of May to detect neutralising antibodies to SARS-CoV-2 as a marker of past infection and epidemic progression. AimTo determine if sera from blood bank donors can be used to track the emergence and progression of the SARS-CoV-2 epidemic. MethodsA pseudotyped SARS-CoV-2 virus microneutralisation assay was used to detect neutralising antibodies to SARS-CoV-2. The study group comprised samples from 3,500 blood donors collected in Scotland between the 17th of March and 19th of May, 2020. Controls were collected from 100 donors in Scotland during 2019. ResultsAll samples collected on the 17th March, 2020 (n=500) were negative in the pseudotyped SARS-CoV-2 virus microneutralisation assay. Neutralising antibodies were detected in 6/500 donors from the 23th-26th of March. The number of samples containing neutralising antibodies did not significantly rise after the 5th-6th April until the end of the study on the 18th of May. We find that infections are concentrated in certain postcodes indicating that outbreaks of infection are extremely localised. In contrast, other areas remain comparatively untouched by the epidemic. ConclusionThese data indicate that sero-surveys of blood banks can serve as a useful tool for tracking the emergence and progression of an epidemic like the current SARS-CoV-2 outbreak.

Estimating false-negative detection rate of SARS-CoV-2 by RT-PCR (preprint)

Reverse transcription-polymerase chain reaction (RT-PCR) assays are used to test patients and key workers for infection with the causative SARS-CoV-2 virus. RT-PCR tests are highly specific and the probability of false positives is low, but false negatives can occur if the sample contains insufficient quantities of the virus to be successfully amplified and detected. The amount of virus in a swab is likely to vary between patients, sample location (nasal, throat or sputum) and through time as infection progresses. Here, we analyse publicly available data from patients who received multiple RT-PCR tests and were identified as SARS-CoV-2 positive at least once. We identify that the probability of a positive test decreases with time after symptom onset, with throat samples less likely to yield a positive result relative to nasal samples. Empirically derived distributions of the time between symptom onset and hospitalisation allowed us to comment on the likely false negative rates in cohorts of patients who present for testing at different clinical stages. We further estimate the expected numbers of false negative tests in a group of tested individuals and show how this is affected by the timing of the tests. Finally, we assessed the robustness of these estimates of false negative rates to the probability of false positive tests. This work has implications both for the identification of infected patients and for the discharge of convalescing patients who are potentially still infectious.

Fundamental principles of epidemic spread highlight the immediate need forlarge-scale serological surveys to assess the stage of the SARS-CoV-2 epidemic (preprint)

The spread of a novel pathogenic infectious agent eliciting protective immunity is typically characterised by three distinct phases: (I) an initial phase of slow accumulation of new infections (often undetectable), (II) a second phase of rapid growth in cases of infection, disease and death, and (III) an eventual slow down of transmission due to the depletion of susceptible individuals, typically leading to the termination of the (first) epidemic wave. Before the implementation of control measures (e.g. social distancing, travel bans, etc) and under the assumption that infection elicits protective immunity, epidemiological theory indicates that the ongoing epidemic of SARS-CoV-2 will conform to this pattern. Here, we calibrate a susceptible-infected-recovered (SIR) model to data on cumulative reported SARS-CoV-2 associated deaths from the United Kingdom (UK) and Italy under the assumption that such deaths are well reported events that occur only in a vulnerable fraction of the population. We focus on model solutions which take into consideration previous estimates of critical epidemiological parameters such as the basic reproduction number (R0), probability of death in the vulnerable fraction of the population, infectious period and time from infection to death, with the intention of exploring the sensitivity of the system to the actual fraction of the population vulnerable to severe disease and death. Our simulations are in agreement with other studies that the current epidemic wave in the UK and Italy in the absence of interventions should have an approximate duration of 2-3 months, with numbers of deaths lagging behind in time relative to overall infections. Importantly, the results we present here suggest the ongoing epidemics in the UK and Italy started at least a month before the first reported death and have already led to the accumulation of significant levels of herd immunity in both countries. There is an inverse relationship between the proportion currently immune and the fraction of the population vulnerable to severe disease. This relationship can be used to determine how many people will require hospitalisation (and possibly die) in the coming weeks if we are able to accurately determine current levels of herd immunity. There is thus an urgent need for investment in technologies such as virus (or viral pseudotype) neutralization assays and other robust assays which provide reliable read-outs of protective immunity, and for the provision of open access to valuable data sources such as blood banks and paired samples of acute and convalescent sera from confirmed cases of SARS-CoV-2 to validate these. Urgent development and assessment of such tests should be followed by rapid implementation at scale to provide real-time data. These data will be critical to the proper assessment of the effects of social distancing and other measures currently being adopted to slow down the case incidence and for informing future policy direction.